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ABSTRACT 



Context. High-mass microquasars are binary systems consisting of a massive star and an accreting compact object from 
which relativistic jets are launched. There is considerable observational evidence that winds of massive stars are clumpy. 
Individual clumps may interact with the jets in high-mass microquasars to produce outbursts of high-energy emission. 
Gamma-ray flares have been detected in some high-mass X-ray binaries, such as Cygnus X-1, and probably in LS 5039 
and LS 1+61 303. 

Aims. We predict the high-energy emission produced by the interaction between a jet and a clump of the stellar wind 
in a high-mass microquasar. 

Methods. Assuming a hydrodynamic scenario for the jet-clump interaction, we calculate the spectral energy distribu- 
tions produced by the dominant non-thermal processes: relativistic bremsstrahlung, synchrotron and inverse Compton 
radiation, for leptons, and for hadrons, proton-proton collisions. 

Results. Signiflcant levels of emission in X-rays (synchrotron), high-energy gamma rays (inverse Compton), and very 
high-energy gamma rays (from the decay of neutral pious) are predicted, with luminosities in the different domains in 
the range ~ erg s The spectral energy distributions vary strongly depending on the speciflc conditions. 

Conclusions. Jet-clump interactions may be detectable at high and very high energies, and provide an explanation for 
the fast TeV variability found in some high-mass X-ray binary systems. Our model can help to infer information about 
the properties of jets and clumpy winds by means of high-sensitivity gamma-ray astronomy. 

Key words. Gamma-rays: theory - X-rays: binaries - Radiation mechanisms: non-thermal 



1. Introduction 

The mass loss in massive and hot stars is understood to 
occur by means of supersonic inhomogeneous winds struc- 
tured as very dense and smaU clumps embedded in large 
regions of tenuous plasma. This idea is supported by con- 
siderable observational evidence of the clumpy structure of 
these winds (e.g. Puis et al. 2006, Owocki & Cohen 2006, 
Moffat 2008). However, as a consequence of the high spa- 
tial resolution necessary for a straightforward detection of 
the clumps, all the evidence is indirect. Hence, properties 
of clumps, such as size, density, and number, arc not well- 
known. 

Some massive stars are accompanied by a compact ob- 
ject, to which they transfer matter. This process leads to 
the formation of an accretion disk around the compact 
object and, in high-mass microquasars (HMMQs), to the 
generation of bipolar relativistic outflows (e.g. Mirabel & 
Rodriguez 1999). 

Non-thermal emission has been observed in microquasar 
jets from radio (e.g. Ribo 2005) to X-rays (e.g. Corbel et al. 
2002). At higher energies, gamma-ray radiation could also 
be produced in jets (e.g. Bosch- Ramon et al. 2006). A TeV 
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flare was detected by MAGIC from the HMMQ Cygnus X- 
1 (Albert et al. 2007). Transient gamma-ray events may 
also have been detected from the high-mass X-ray bina- 
ries LS 5039 and LS 1+61 303 by HESS (Aharonian et al. 
2005) and MAGIC (Albert et al. 2006), respectively, as sug- 
gested by Paredes (2008). In addition, the HMMQ Cygnus 
X-3 might have been observed flaring in the GeV range 
by AGILE (ATels 1492, 1547 and 1585, see however Atel 
1850). Both this instrument and Fermi also found several 
transient GeV sources in the Galactic plane without known 
counterpart (ATel 1394, Abdo et al. 2009). AU this strongly 
variable gamma-ray emission may have a similar origin. For 
instance, the interaction between the jet and the stellar 
wind of the companion star could produce TeV flares (e.g., 
Romero et al. 2003, Romero & Orellana 2005, Romero et 
al. 2007, Albert et al. 2007, Perucho & Bosch-Ramon 2008, 
Owocki et al. 2009). 

In this paper, we propose a model to explain these 
gamma-ray flares, based on the interaction between the 
jets of a HMMQ with wind inhomogeneities. The clumps 
can eventually penetrate the jet, leading to transient non- 
thermal activity that may release a significant fraction of 
the jet kinetic luminosity in the form of synchrotron, inverse 
Compton (IC), and proton- proton (pp) 7r°-decay emission. 

This work is organized as follows. In the next section, 
we describe the main characteristics of the scenario adopted 
and estimate the relevant timescales of the jet-clump inter- 
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Fig. 1. Sketch of a HMMQ with clumpy stellar wind, 
adapted from Romero et al. (2007). 

action; in Sect. 3, we study the acceleration of particles and 
estimate the non-thermal radiative losses; in Sect. 4, we 
describe the calculation of the non-thermal emission and 
present the main results; and finally, in Sect. 5, we draw 
some conclusions. CGS units are used consistently through- 
out the paper. 

2. The physical scenario 

To study the interaction between a clump of the stellar wind 
and a jet in a HMMQ, we adopt a scenario with similar 
characteristics to the binary system Cygnus X-1. We fix 
the separation between the compact object and the massive 
star to be a = 3 X 10^^ cm (0.2 UA). For the star luminosity 
and temperature, we assume that — 10^^ erg s~^ and 
T^, = 3 X 10* K, and the stellar mass loss rate is adopted to 
be Mi, = 3 X 1O~^M0 yr~^, with a terminal wind velocity 
Ww ^ 2.5 X 10* cm s~^. A sketch of the scenario is presented 
in Fig. [B 

2.1. Clump model 

The clump is assumed to be spherical and homogeneous. 
Given the uncertainties in the clump parameters, we con- 
sider two values of its size: i?c — 10^° and 10^^ cm (i.e., 
'--^ 3 X (10~^ — 10~^)a). For effective jet penetration, a 
large density contrast between the clump and the jet is 
required. Therefore, we assume dense clumps with ric = 
10^^ cm"'', which correspond to a porous clumpy wind with 
a filling factor (clump versus interclump volume ratio) of 
/ = Mi,/ ^'Ko^mpVc'ric — 0.005, where rup is the mass of the 
protorQ. We assume that the velocity of the clumps equals 
the velocity of the wind, i.e., Vc = 2.5 x 10* cm s~^. The 
temperature of the clumps is taken as Tc = 10'' K (Krticka 

^ See Owocki & Cohen (2006) for the concept of a porous 
wind. 



Table 1. Adopted parameters in this work. 



Parameter description [unitsj 


Clump 


Jet 


Radius [cm] 


10^" - 10^^ 


1.5 X 10'' 


Velocity [cm s^""^] 


2.5 x 10* 


10^° 


Density [cm^'^] 


1012 


4.7 X 10^ 


Binary system 


System size [cm] 


3 X 10'^ 




Star luminosity [erg s^^] 


^039 




Star temperature [K] 


3 X 10* 




Mass loss rate [M© yr~^] 


3 X lO"** 




Wind velocity [cm s~^] 


2.5 X 10** 





& Kubat 2001), which is moderately lower than the tem- 
perature at the surface of the massive star. 

2.2. Jet model 

Radio observations demonstrated that the jets of MQs are 
strongly coUimated (Miller- Jones, Fender & Nakar 2006). 
We assume that the jet radius is one tenth of the jet 
height, i.e. Rj{z) — 0.1 z, which corresponds to a semi- 
opening angle oi (f) = 6°. The expansion velocity results 
in Vexp = 0.1 Vj, where Vj is the velocity of the jet. This 
expansion velocity for a free-expanding hidrodynamical jet 
implies a jet base Mach number Vj/v^xp ^ 10, i.e. a 
strongly supersonic outflow. We consider a jet dynamically 
dominated by cold protons with a weakly relativistic bulk 
velocity Vj = 0.3 c, i.e. a Lorentz factor 7j = 1.06. We 
neglect the curvature of the jet produced by its interac- 
tion with the stellar wind. In HMMQ jets of luminosity 
> 10^^ erg s"-', the geometry should not be strongly mod- 
ified (Perucho & Bosch- Ramon 2008), although this effect 
could be important in systems such as Herbig-Haro jets 
interacting with a stellar wind, as studied by Raga et al. 
(2009). 

The kinetic luminosity of the jet is taken to be Lkin ~ 
3 X lO'^^ erg s^^, similar to that of Cygnus X-1 (e.g., Gallo 
et al. 2005, Russell et al. 2007). Using the equation 

-^kin = o"j(7j - l)mp rij Vj, (1) 

where ctj — ttR^ is the cross-section of the jet, we can es- 
timate the density rij of the jet material in the laboratory 
reference frame (LRF). At the height of the jet-clump in- 
teraction, which is taken to occur at z = a/2, we obtain 
Tij = 4.7 X 10^ cm"''. Thus, the ratio of the clump to the jet 
densities is x = 2.1 x 10*. This parameter will be very rele- 
vant to the jet-clump interaction estimates. The parameter 
values of the clump and the jet are shown in Table [TJ We 
assume an ideal gas equation of state, P = nkT, for both 
the clump and the jet, where n is the gas density. 

2.3. Dynamics of the interaction 

Clumps are formed in the wind acceleration region within 
one stellar radius of the surface of the star (e.g. Puis et al. 
2006) and some of them reach the jet/wind interface. The 
very large inertia of the clumps, linked to the exceptionally 
high density contrast x, allows them to cross the boundary 
of the jet and fully penetrate into it. A large value of x 
is also required to ensure that the clump is not strongly 
affected when penetrating into the jet. 
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To study the physical processes of the interaction, we 
consider the collision of a single clump with the jet. For 
simplicity, we assume that, on the relevant spatial scales 
of the interaction, the jet is cylindrical. In the context of 
this work, thermal conduction, clump expansion, magnetic 
fields and gravitational forces are not dynamically relevant 
and will be neglected. 

In the LRF, the clump will take a time tc, or jet pene- 
tration time, to fully enter the jet. This provides the first 
relevant timescale 



tr. 



2Rc/vc, 



(2) 



which is tc ~ 80 and ~ 800 s, for = lO^"^ and 10" cm, 
respectively. In addition, the clump crosses the jet roughly 
at the wind velocity in the jet-crossing time 



2 i?i 



= 1.2 X 10-^ s. 



(3) 



From the moment when the clump interacts with the jet, 
the ram pressure exerted by the latter produces a shock in 
the clump, which propagates in the direction of the jet mo- 
tion. Assuming that a significant fraction of the momentum 
flux is transferred to the clump, the clump-crossing time, 
which is the characteristic timescale of the jet-clump inter- 
action, can be defined to be 



tr. 



2 Rr 2 i?,VX 



where 



J2_ 



(4) 



(5) 



is the velocity of the shock moving through the clump, 
derived by equating the jet/clump ram pressures. The 
timescale tec is - 3 x 10^ and 3 x 10^ s for i?c = 10^° 
and 10^^ cm, respectively. 

A shock (the bow shock) is also formed in the jet when 
its material collides with the clump. We assume that the 
bow-shock region width (or clump/bow-shock separation 
distance) is x 0.2i?c (van Dyke & Gordon 1959), and 
thus the time required to reach the steady state regime is 



ths 



0.2Rc 



(6) 



where wj ps ^ Vj/A is the jet postshock velocity in the LRF. 
Therefore, tbs ^ (5/2)icc/Vx, i.e. ibs < ^cc- In Fig. H a 
sketch of the situation is shown. 

Once the clump is inside the jet, the latter transfers 
momentum to the clump, accelerating it to the background 
velocity, Vj . The acceleration can be obtained from the force 
exerted by the jet ram pressure and the clump mass (e.g.. 
Fragile et al. 2004): 



xRc 



(7) 



The jet accelerates the clump to the velocity Vj in a time 



t^ 



9 



(8) 



which is much longer than t^c . When one fluid exerts a force 
against another fluid of different density, the hydrodynam- 
ical Rayleigh- Taylor (RT) instability eventually develops. 




Bow-shock 



Jet 



Clump 




Massive 
Star 



Compact Object 
Fig. 2. Sketch of the jet-clump interaction. 



leading to the perturbation and potential disruption of the 
clump. The timescale for this instability is 



tRT ^ v^/g 



(9) 



where I is the instability length in the perturbed region. 

After the bow shock is formed, the jet material sur- 
rounds the clump and we have two fluids that have a large 
relative velocity. This situation leads to Kelvin-Helmholtz 
(KH) instabilities. The timescale for a KH instability to 
develop is 



^KH 



Vrol 



(10) 



where Vi-d is the relative velocity between the jet material 
and the clump material (we assume Vrd ~ Vj). In addition, 
taking into account that we are interested in instability 
lengths similar to the clump size, i.e. I Rc, we find that 

^RT ^ tKH ~ tec- 

According to the timescales estimated in the previous 
paragraphs, the clump can fully enter into the jet if tc ~ 
tcc/5. At this stage, we do not consider the penetration 
of the clump into the jet, and assume that the former is 
completely inside the latter (i.e., the system has cylindrical 
symmetry) . The bow shock is formed in a time much shorter 
than tc and ice (^bs ~ tcc/y/x)- We note that the clump 
might not escape the jet if tj > tRT/KH- In such a case, the 
clump will be destroyed inside the jet. However, numerical 
simulations show that the instability timescales are longer 
than the clump crossing time by a factor of a few (e.g., 
Klein, McKee & Colella 1994), i.e., tRx/KH > ice- 
Regarding the shock properties and given the particu- 
lar characteristics of our scenario, the shock in the clump 
is strong, radiative, and slow, whereas the bow shock is 
strong as well, but adiabatic and fast. For these reasons, 
the shocked and heated material of the clump will radiate 
a non-negligible part of the energy tranferred by the jet. 
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2.3.1. Thermal emission from the clump 

To estimate the density, nps(C), and temperature, Tps(C), 
of the shocked clump at a distance C from the shock, we as- 
sume pressure equilibrium between the clump and the jet 
shocked plasma and use the equation of energy flux conser- 
vation, taking into account the main channels of thermal 
radiative losses: A(T) = 7 x 10"^^ T-°-^ erg cm^^ g-i. The 
timescale of thermal cooling equals 



ith = 3 X 10' 



1.6 



^ps(C) 



(11) 



The temperature and density in the adiabatic zone of the 
post-shock region are 8.5 x 10^ K and 4 x 10^^ cm~^, re- 
spectively. From these values, the shocked clump material 
cools in a time ^th ~ 10 s, which is shorter than tec, which 
means a distance Cth ~ ithWcs/4 ~ 2 x 10^ cm from the 
shock front. Since ^th < 2i?c, the shock in the clump is 
radiative. At distances greater than ^th, the density of the 
shocked matter can reach values as high as 10^^ cm~^. 

Although we are mainly interested in the high-energy 
emission produced by the jet-clump interaction, we esti- 
mate for comparison the free-free emission generated by the 
shocked and heated clump. Considering rps(C) and nps(C), 
we estimate the free- free luminosity by integrating the emis- 
sivity (Lang 1999) along the shocked clump, as demon- 
strated in Zhekov & Palla (2007). The corresponding ther- 
mal luminosities are Lth ~ 5 x 10'^° and 5 x 10^^ erg s~^ 
for Rc = 10^° and 10^^ cm, respectively, peaking in the 
soft X-rays. In Fig. [3 the thermal radiation of the clump 
is shown together with the non-thermal emission. 

In contrast to what occurs in the clump, the bow shock 
is adiabatic and fast. For this reason, it is a propitious place 
to accelerate particles to relativistic energies. 



3. Particle acceleration and non-thermal cooling 

3.1. Particle acceleration 

In the presence of a shock and magnetic field, non- 
relativistic diffusive (Fermi I) shock acceleration can occur 
(e.g. Drury 1983). In the linear limit of this theory and 
considering Bohm diffusion, electrons and protons will be 
accelerated to an energy E^^p (where e stands for electrons 
and p for protons) in a time 



^acc — ^ 



E, 



qBc' 



(12) 



where 7] ~ (8/3)(c/t>s)^ for perpendicular shocks (Protheroe 
1999). In Eq. p^ . Vg is the shock velocity (taken to be equal 
to Vj), B is the magnetic field in the acceleration region, c 
is the speed of light, and q is the electron charge. Thus, to 
obtain efficient acceleration of particles, i.e., a short taccj a 
high value of B and a strong and fast shock are necessary. 

To study the jet-clump interaction, we consider two val- 
ues of the magnetic field in the jet shocked (bow-shock) re- 
gion, B\js- First, we consider a value of the magnetic field 
obtained by assuming that the magnetic energy density is 
10 % of the plasma internal energy density downstream. 
We adopt a sub-equipartition value to ensure that Bbs is 
dynamically negligible with respect to matter. Fixing 



87r 



O.Iud 



(13) 



(14) 



where Up is the energy density of the shocked jet gas, given 
by 

3„ 9 

we obtain B^s 150 G. This magnetic field, plus Vg — Vj, 
yield iacc ~ lO^^Se.p s. 

In addition, we adopt Sbs = 1 G to check the impact 
on our results of a magnetic field significantly weaker than 
in the sub-equipartition case. 

3.2. Non-thermal cooling processes 

Relativistic leptons lose their energy by different non- 
thermal radiative processes, such as synchrotron radiation, 
IC scattering and relativistic bremsstrahlung. On the other 
hand, relativistic protons can also lose energy by means of 
pp interactions. Finally, the shocked plasma can produce 
thermal radiation if the density is high enough, as shown 
in Sect. \^XT\ 

Relativistic electrons that lose their energy by syn- 
chrotron radiation have the following cooling time (e.g. 
Ginzburg & Syrovatskii, 1964) 



4.1 X 10^ 

B^E, 



(15) 



Synchrotron radiation is the most efficient radiative process 
in the bow-shock region in the case of = 150 G, where 
Wnch-2x 10-2i;-i s. 

At the interaction height considered in this work, Zjnt = 
a/2 = 1.5 X 10^^ cm, the energy density of the photons from 
the star is Uph = 2.4 x 10^ erg cm^'^, the typical photon 
energy being eo ^ 10 eV. For y = eoi?e/(5.1 x 10^ eV)^ > 1, 
i.e. Ee > 2.6 X 10^'^ eV, the IC interaction occurs in the 
Klein-Nishina (KN) regime. A formula for the IC cooling 
time valid in both a Thompson (Th) and KN regime in a 
photon field with a narrow energy distribution is e.g. Bosch- 
Ramon & Khangulyan (2009) 

^ _ 6.1xlOi^eo (l + 8.3j/) (1 + 1.3y^) 

Mph ln(l + 0.2y) (l + 0.5y-f 1.3?/) ^ > 

For i? ^ 1 G, IC is the dominant loss leptonic channel. 

According to the low particle density of the shocked jet, 
nbs = 4nj = 2 X 10^ cm~'^, relativistic bremsstrahlung 
losses are negligible in the bow-shock region, the cool- 
ing time for a completely ionized medium being (e.g. 
Blumenthal & Gould, 1970): 



1.4 X 10 



16 



(i^(;;!^-)+o.36) 



(17) 



where Z is the atomic number, and in our case Z = 1. 

Regarding hadronic emission, 7-rays are produced if rel- 
ativistic protons interact with nuclei through inelastic col- 
lisions. As in the case of relativistic bremsstrahlung, the 
proton cooling time due to pp interactions depends on the 
density n (e.g. Aharonian & Atoyan, 1996) related 



tpp 



2 X 10^ 



(18) 



In the bow-shock region, this process is negligible. 
Otherwise, if relativistic protons accelerated in the bow 
shock penetrate into the clump, pp interactions can become 
an efficient process to generate gamma-rays. 
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log(E, [eV]) 

Fig. 3. Acceleration and radiative loss (synchrotron and 
IC) time for electrons in the bow-shock region. The advec- 
tion and diffusion times are shown for Rc = 10^° (thin line) 
and 10^^ cm (thick line). This figure corresponds to the case 
Bbs = 1 G. 
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Fig. 4. Acceleration and radiative loss (synchrotron and 
IC) time for electrons in the bow-shock region. The advec- 
tion and diffusion times are shown for Rc = 10^'' (thin line) 
and 10^"'^ cm (thick line). This figure corresponds to the case 
Bbs = 150 G. 



3.3. Maximum energies 

Taking into account energy gains and losses, the maxi- 
mum energy achieved by particles accelerated in the bow 
shock can be easily estimated. Concerning energy losses, 
we consider radiative cooling (described above) and the 
escape of particles. The latter takes into account the ad- 
vection of relativistic particles by downstream bow-shock 
material (tadv ~ Rc/vips) and the diffusion of particles 
(idiff ^ x'^/2Db, where x = 0.2 R^ and Db = -Ee,pc/3gi?bs 
is the diffusion coefficient, which is assumed to be the Bohm 
one). The corresponding timescale is 

Tcsc = min{tdiff , tadv} • (19) 

Given the small thickness of the bow-shock region, we as- 
sume that particles that escape via diffusion go to the 
clump. 

For electrons, the maximum energy is constrained by the 
escape of particles via diffusion (Bbs = 1 G) and by syn- 
chrotron radiation (Bbs = 150 G), as it is shown in Figs. 
[3] and m From these figures, we determine that the most 
relevant radiative process in the bow-shock region can be 
IC scattering if Bbs = 1 G, and synchrotron radiation if 
Bhs = 150 G. On the other hand, the maximum energy for 
protons accelerated in the bow shock is constrained by the 
Hillas criterion (Hillas 1984), i.e. when the proton gyrora- 
dius becomes equal to the size of the acceleration region: 

E^^^ = xq Bbs = 0-2 RcqB^,. (20) 

In Table [51 maximum energies for electrons and protons of 
different values of Bbs and Rc are shown. 

4. Production of gamma-rays and lower energy 
radiation 

We calculated the spectral energy distribution (SED) of 
the emission produced by the most relevant non-thermal 



Table 2. Maximum energies achieved by acelerated par- 
ticles in the bow shock. The values are in eV units and 
represent the different values of Rc and B studied in this 
work. 



-Rc [cm] 


10^" 


10^" 


10" 


10" 


Bbs [G] 


1 


150 


1 


150 


B^"" [eV] 


1.5 X 10"" 


8 X 10'' 


1.5 X 10'^ 


8 X 10" 


-B^"" [eV] 


6 X 10" 


9 X 10^^ 


6 X 10^^ 


9 X 10^'' 



radiative processes. In the bow-shock region, we considered 
synchrotron and IC radiation, and in the clump, we also 
considered relativistic bremsstrahlung and pp. 

4.1. Distribution of relativistic particles 

We assume an injected population of relativistic particles 
(electrons and protons) in the bow-shock region that follows 
a power-law energy distribution of the form 

Qe,p{Ee,p) = i^e,p e^" E^^P / Kp"") ■ (21) 

The index F is fixed to 2, typical of linear diffusive shock 
acceleration, and we add an exponential cut-off at high en- 
ergies. The normalization constant K^^p is determined by 
assuming that 25 % of the jet luminosity incident on the 
bow shock, i.e. ~ ((7c/CTj)Lkin (where Cc = t^Rc is the clump 
effective cross section p) , is converted into power of the ac- 
celerated particles. We note that the predicted fluxes scale 
linearly with the adopted non-thermal fraction. 

To estimate the particle energy distribution N^^p, we 
solve the kinetic equation in the one-zone model approxima- 



^ We neglect here the region where the bow shock be- 
comes strongly oblicuous, and focus on where this shock is the 
strongest, i.e. right in front of the clump. 
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tion for the bow-shock region (e.g. Ginzburg 
1964) 



d 



dt 



N, 



P ' 



Syrovatskii 



(22) 



where t is the time and Pe.p — —dE^.p/dt is the energy loss 
rate for particles. 

The relativistic leptons reach the steady state well be- 
fore the shock has crossed the clump. As shown in Figs. [3] 
andUl the most energetic electrons can diffuse to the clump 
(i?bs = 1 G) or lose their energy inside the bow-shock region 
by synchrotron and IC radiation {B^s = 150 G). However, 
particles downstream with low energies escape advected in 
the shocked material of the jet before cooling radiatively, 
producing a break in the energy spectrum of particles. By 
equating the advection and synchrotron loss times, we can 
estimate the break energy in the case with Bbs — 150 G, 
which results in = 3 x 10^ and £'b = 3 x 10*^ eV, 
for i?c = 10^° and 10^^ cm, respectively. For the case of 
Bbs = 1 G, the break energy is determined by the advec- 
tion and diffusion times, given E\^ — 0.3E™^^. The electrons 
with energies E,, > E^ can reach the clump and radiate in- 
side it. The energy distribution of these electrons is 



(23) 



where Ql is the part of Qe that corresponds to 0.3 E'^^^ < 
E, < i?--, and ieooi = l/it-i + t-i + t-^^ + trl). The 
ionization cooling term, tion ~ 2 x 10^^(£'e/'^) s, comes 
from the high densities in the clump. 

As noted in Sect. 13. 2i relativistic protons do not suffer 
significant pp losses in the bow-shock region {tditf ^ tpp)- 
These protons can also reach the clump if they are not 
advected by the shocked material of the jet. By assuming 
that idiff < iadv7 the minimum energy necessary to reach 
the clump is 

Ef'' = 0.025£;^"'^'^ , (24) 

and the maximum energies of these protons are determined 
by the Hillas criterion and shown in Table [2l On the other 
hand, to confine relativistic protons in the clump, the mag- 
netic field must be > 10'^ G, much stronger than the ex- 
pected clump magnetic fielqj, i.e., the protons will cross 
the entire clump in a time ~ Rc/c < tpp. We note that Rc 
should be the shocked clump width if the clump has already 
been affected by the shock. The distribution of relativistic 
protons in the clump is 



Np{Ep) 



Rc 



Qp{Epi 



(25) 



With the steady distributions of relativistic electrons in the 
bow shock, Ne{Ee), and of both electrons and protons in the 
clump, N^{Ee) and Np{Ep), respectively, we can calculate 
the SEDs of the radiation produced by these particles. 

4.2. Non-thermal emission from the bow shock 

Since the energy density of the synchrotron emission is 
lower than that of the magnetic and stellar fields, syn- 
chrotron self-Compton processes will be neglected in our 
calculations. 



^ Assuming equipartition, magnetic fields ~ 100 G would be 
expected. 



Synchrotron radiation is computed using the standard 
formulae given in Blumenthal & Gould (1970). Assuming 
that the distribution of relativistic electrons in the bow- 
shock region is isotropic and moves with a non-relativistic 
advection speed, we calculate the synchrotron SED as 

eisynch(e) = e / ' N,{E,) P{Ee,e) dE, , (26) 

where 

P{E,, e) = 1034i?b.(e/£;c)('/'^ exp(-e/i?c) s'^ (27) 

is the power function. The magnitude Ec{Ee) = 5.1 x 
10~^ B El erg is the characteristic energy of the produced 
photons. 

For a thermal distribution of target photons, we can 
estimate the IC emission from 

,Lic(6) = Uphce Ar^(£;^)^^^(^I£l!^di?e, (28) 

J fimin atje 

where da{xic,iQ, E^) / dEf. is a parametrized differential 
cross section for an isotropic target photon field, valid in 
both Th and KN regimes (Blumenthal & Gould 1970). This 
magnitude depends on the adimensional parameter 



(a-) (i-t) 



(29) 



for m1c*/{AE'^) < xic < 1, where eo is the energy of the 
target photons. 

We calculated eLsynch(e) and eLic(e) for different values 
of i?bs and Rc, and the results are presented in Figs. [5] 
and [51 As seen in these figures, the synchrotron component 
is more luminous than the IC one in the cases of Bbs = 
150 G, reaching bolometric luminosities of igynch ~ 10'^^ 
and 2x 10'^^ erg s~^ for Rc = 10^° and 10^^ cm, respectively. 
In contrast, for By,s = 1 G, the dominant radiative process 
is IC scattering, which leads to bolometric luminosities of 
Lie ~ 2 X 10^^ erg s'^ and ~ lO^^ erg s'^ for Rc = 10^° 
and 10^^ cm, respectively. The maximum energies achieved 
by photons can be as high as E'^\^^ ^ 1 TeV. We note the 
spectral break in the synchrotron and IC emission for the 
case i?bs — 150 G. This feature is produced by the advection 
of particles out from the bow-shock region. In addition, the 
impact of KN losses hardening the electron spectrum can 
also be seen in the spectral shape of the synchrotron and 
IC emission for Rc — 10^^ cm and i?bs = 1 G. 

We accounted for 7 — 7 absorption in the stellar field 
of the massive star. Since we did not focus on a particular 
binary/observer system geometry, we assumed an isotropic 
target photon field (and hence also neglected angular ef- 
fects in the IC calculations) . As seen in Figs. [5] and [6l the 
attenuation in flux can be of several orders of magnitude at 
energies of hundreds of GeV. In only some cases of specific 
geometries in the 7 — 7 interaction, can the attenuation 
be very much reduced (see e.g. Khangulyan, Aharonian & 
Bosch- Ramon 2008). 

In this work, we have focused on the high-energy emis- 
sion. Radio emission is also produced by electrons that leave 
the bow shock well before they lose most of their energy. 
These particles may radiate in the radio band further down 
the jet and their treatment is beyond the context of this 
study. Because of this, we have not considered the effect of 
synchrotron self-absorption in the radio spectrum, and we 
do not make predictions in this energy range. 
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Fig. 5. Computed synchrotron (dotted line) and IC (dashed 
hne) SEDs for the emission produced in the bow-shock re- 
gion for the case of i?c — 10^° cm with Bbs — 1 (thin hne) 
and 150 G (thick hne). The production and the 7 — 7 ab- 
sorbed curves at very high energies are shown. 
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Fig. 6. The same as in Fig.[5l but for the case i?c = 10^^ cm. 



4.3. Non-thermal emission from the clump 

The most energetic particles accelerated in the bow-shock 
region can penetrate the clump. In the case of protons, 
uncooled particles with Ep > 0.025E™'^^ can diffuse up to 
the clump and radiate only a part of their energy from 
there. On the other hand, for Bbs = 1 G (< Bequipartition), 
electrons with > 0.3i?™^'^ will radiate their energy in 
the clump. 

The relativistic protons that diffuse up to the clump col- 
lide with cold protons there, generating both neutral (tt*^) 
and charged (tt^) pions that decay to 7-rays and leptons, 
respectively. 

The emissivity of tt^ is given by 



{Ep)Np{Ep)dEp , 



where app{Ep) is the cross section of pp interactions es- 
timated by Kelner et al. (2006). With knowledge of the 
emissivity of pions, the specific luminosity of the photons 
produced by Tr^'-decay is 



-.dE^ 



(31) 



(30) 



where i;^'" = e -I- mlc'^/i'ie) and i;;"^'^ - E^'''''/6. The 
estimated emission is shown in Fig. [71 along with the con- 
tributions of the other processes. The emission reaches lu- 
minosities of as high as Lpp ~ 10^^ erg s~-^ {Rc — 10^^ cm; 
Fig. [7]). For the adopted values of _Bbs and E^'^^, the 7 — 7 
absorption reduces the final fluxes and can lead to very 
different spectral shapes. 

Secondary pairs with an effective low energy cut-off at 
^ 0.1 £'p"" will be injected inside the clump by the de- 
cay of TT^. These pairs radiate most of their energy inside 
the clump by synchrotron, IC scattering, and relativistic 
bremsstrahlung. In general, primary electron emission from 
both the bow-shock region and the clump (see the next 
paragraph) is more significant than that of this secondary 
component (see Bosch- Ramon, Aharonian & Paredes 2005, 
and Orellana et al. 2007 for discussion of secondary pairs 
in the context of clouds and jets, respectively). Finally, 
we note that very high-energy neutrinos with luminosities 
~ Lpp would also be generated (e.g., Aharonian et al. 2006; 
Reynoso & Romero 2009). 

To estimate the radiation produced by electrons accel- 
erated in the bow-shock, we assume two values of the mag- 
netic field of the clump: Be ^ 1 and 100 G, the latter being 
similar to the equipartition clump magnetic field. 

The synchrotron and the IC radiation were calculated 
using the same equations given in Sect. 14.21 Regarding 
the relativistic bremsstrahlung, we used the standard for- 
mulae given in Blumenthal & Gould (1970). As shown in 
Fig. [8l synchrotron radiation is higher than IC in the case 
of i?c = 100 G, reaching a luminosity ~ 10'^'* erg s~^ 
{Rc = 10" cm) at e ~ 0.1 MeV. For Be = 1 G, on the 
other hand, the absorbed IC component reaches a similar 
luminosity to that of the synchrotron one, ~ 10^^ erg s~^. 
Relativistic bremsstrahlung is negligible in all cases. The 
SEDs for Rc = 10^° cm are similar to those shown in Fig. 
[51 but the luminosities are oc and therefore two orders 
of magnitude less. We note that the non-thermal emission 
from the clump is similar to, but of slightly lower intensity 
than the bow shock one, being the spectra harder. We note 
that the radio emission produced by leptons in the clump is 
strongly suppressed by ionization losses and the high low- 
energy cut-off of the injected electrons. For clarity, we have 
not plotted the leptonic emission from the clump together 
with the other components, but they can be compared us- 
ing Figs [3 and [51 



5. Discussion and summary 

We explored the main physical processes and the nature 
of the radiation produced by the interaction between the 
jet of a HMMQ and a clump of the wind of the compan- 
ion star. The penetration of a clump in the jet produces 
two shocks. One shock quickly reaches a steady state in 
the jet, forming a bow-like shock. The other shock prop- 
agates through the clump driven by pressure equilibrium 
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Fig. 7. SEDs for different values of Bbs and i?c; the curves of both absorbed and unabsorbed (thin lines) IC and -p-p 
radiation are shown. 



in the jet/clump contact discontinuity. The bow shock is 
adiabatic and fast, and particles can be accelerated to very 
high energies, leptons efficiently emitting synchrotron and 
IC radiation. Otherwise, the shock in the clump is slow 
and radiative. It is a non-efficient accelerator but a ther- 
mal emitter as well as a good target for proton-proton col- 
lisions. If electrons from the bow-shock region enter the 
clump, they could also radiate energy efficiently there via 
synchrotron and IC processes. The SEDs of the mentioned 
radiation components have been computed in the context of 
a microquasar with parameters similar to those of Cygnus 
X-1, showing that in some cases the interaction between 
one clump and the jet can yield significant amounts of ra- 
diation. The entire SEDs for different cases of clump size 
and magnetic field are shown in Figs. [7] and [51 

As noted above, given the limitations of the one-zone 
approximation, the synchrotron emission has been com- 
puted in the optically thin case, neglecting the impact of 
synchrotron self-absorption. Nevertheless, we note that the 
radio fluxes can be quite high, despite probably being self- 
absorbed right at the bow-shock region. Radio emission 



produced further down the jet by electrons accelerated in 
the bow shock may contribute significantly to the radio 
band, i.e., jet-clump interactions should have associated ra- 
dio fiares. 

At X-rays, the emission is produced by synchrotron 
(bow shock and clump) and thermal radiation (clump). In 
the cases with i?bs ~ 1 G, the thermal emission, reaching 
luminosities of Lth ~ 10'^^ erg s~^, is of higher intensity 
than that of synchrotron from the bow-shock region, but 
not from the clump if = 10^^ cm. Synchrotron emission 
from the bow-shock region dominates in the X-rays band 
for i?bs = 150 G and £synch ^ 10'^^ erg s^^. In a source 
such as Cygnus X-1, these levels of X-ray emission would 
be overcome by the accretion disk radiation. Nevertheless, 
in the case of fainter jet sources at X-rays, such as LS 5039 
and LS I -1-61 303 (Bosch-Ramon et al. 2007, Paredes et 
al. 2007), the synchrotron X-rays produced during the jet- 
clump interaction may be detectable, and even the thermal 
component may be detectable on the top of a non-thermal 
component in certain conditions (large clumps of relatively 
low densities). 
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(and probably also at radio frequencies) of transient activ- 
ity in HMMQs, which are a new tool for studying the winds 
of massive stars. 

Depending on the wind filling factor (or clump density) 
and the clump size/number -TVc- (Owocki & Cohen 2006), 
the flares produced by interactions of clumps with jets in 
HMMQs could be a sporadic phenomenon, for a low num- 
ber of clumps, or may appear as a modulated steady activ- 
ity, for a high number of clumps (see Owocki et al. 2009) . 
For the parameters adopted here, the former case corre- 
sponds to i?c — 10^^ cm and the latter to i?c = 10^° cm. 
Nevertheless, we note that the jet may be disrupted in those 
cases when too many clumps are simultaneously present 
inside the jet. Assuming that jet disruption takes place for 
(Tj < NcXac, for the wind and jet properties adopted in this 
work, the jet could be destroyed if Rc ^ 10^° cm. However, 
more detailed calculations of the dynamics of the jet-clump 
interaction are required to clarify this issue. 



Fig. 8. Non-thermal leptonic emission from the clump, for 
the case Rc — 10^^ cm and for the two assumed values of 
Be, 1 (bottom) and 100 G (top), is shown. The curves for 
both absorbed (thick lines) and unabsorbed (thin lines) IC 
radiation and relativistic bremsstrahlung are shown. 



Inverse Compton scattering in the bow-shock region 
and the clump produces 7-rays up to very high energies 
that dominate the radiation output in cases of relatively 
weak magnetic fields (e.g., i?bs = 1 G). In our calculations, 
the highest luminosity achieved is Lie ~ 10'^^ erg s~^ for 
Rc — 10^^ cm, although 7 — 7 absorption can reduce the 
emission substantially above 100 GeV. Proton-proton colli- 
sions in the clump can also produce 7-rays at energies that 
may be as high as ^ 10^"^ eV (Bbs = 150 G). The maxi- 
mum luminosity obtained by pp is nevertheless quite mod- 
est, Lpp ^ 10"^^ erg s"^ for Rc — 10^^ cm, although denser 
and/or bigger clumps, and more powerful jets may yield 
detectable amounts of photons outside the 7 — 7 absorp- 
tion range (0.1-10 TeV). We recall that specific geometries 
of the binary /observer system plus a high-energy emitter 
far from the compact object may render the attenuation of 
7-rays much smaller (e.g., Khangulyan et al. 2008). 

As a consequence of the characteristics of the interac- 
tion, the expected emission is transient. The flare duration 
is related to the permanence of the clump inside the jet 
dependent on RT and KH instabilities, which can destroy 
the clump. Since the clump is accelerated inside the jet, it 
may be disrupted after several dynamical timescales. If the 
clump were not destroyed, it would eventually leave the jet 
as a thin and moderately hot slab of plasma, since the shock 
in the clump is radiative. Given the typical dynamical and 
crossing timescales, the entire event may have a duration 
of between a few minutes and hours. 

The flares of jet-clump interactions would have associ- 
ated lower (synchrotron, thermal emission) and higher en- 
ergy components (IC, pp emission), which should not be 
correlated with the accretion disk activity. The total level of 
emission, the importance of the different components, and 
the duration of the flares, can provide information about 
the jet power and size, clump size and density, and mag- 
netic fields in the interaction regions (Romero et al. 2007). 
Therefore, besides the jet itself, the clump properties can 
be probed by observations at high and very high energies 
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ABSTRACT 

Context. High-mass microquasars are binary systems consisting of a massive star and an accreting compact object from 
which rclativistic jets arc launched. There is considerable observational evidence that winds of massive stars arc clumpy. 
Individual clumps may interact with the jets in high-mass microquasars to produce outbursts of higli-oucrgy emission. 
Gamma-ray flares have been detected in some high-mass X-ray binaries, such as Cygnus X-1, and probably in LS 5039 
and LS 1+61 303. 

Aims. We predict the high-energy emission produced by the interaction between a jet and a clump of the stellar wind 
in a high-mass microquasar. 

Methods. Assuming a hydrodynamic scenario for the jet-clump interaction, we calculate the spectral energy distribu- 
tions produced by the dominant non-thermal processes: relativistic bremsstrahlung, synchrotron and inverse Compton 

radiation, for loptons, and for hadrons, proton-proton collisions. 

Results. Significant levels of emission in X-rays (synchrotron), high-energy gamma rays (inverse Compton), and very 
high-energy gamma rays (from the decay of neutral pions) are predicted, with luminosities in the different domains in 
the range ~ lO'^^-lO^^ erg s^^. The spectral energy distributions vary strongly depending on the specific conditions. 
Conclusions. Jet-clump interactions may be detectable at high and very high energies, and provide an explanation for 
the fast TeV variability found in some high-mass X-ray binary systems. Our model can help to infer information about 
the properties of jets and clumpy winds by means of high-sensitivity gamma-ray astronomy. 

Key words. Gamma-rays: theory - X-rays: binaries - Radiation mechanisms: non-thermal 



1. Introduction 

The mass loss in massive and hot stars is understood to 
occur by means of supersonic inhomogeneous winds struc- 
tured as very dense and smah clumps embedded in large 
regions of tenuous plasma. This idea is supported by con- 
siderable observational evidence of the clumpy structure of 
these winds (e.g. Puis et al. 2006, Owocki & Cohen 2006, 
Moffat 2008). However, as a consequence of the high spa- 
tial resolution necessary for a straightforward detection of 
the clumps, all the evidence is indirect. Hence, properties 
of clumps, such as size, density, and number, are not well- 
known. 

Some massive stars are accompanied by a compact ob- 
ject, to which they transfer matter. This process leads to 
the formation of an accretion disk around the compact 
object and, in high-mass microquasars (HMMQs), to the 
generation of bipolar relativistic outflows (e.g. Mirabel & 
Rodriguez 1999). 

Non-thermal emission has been observed in microquasar 
jets from radio (e.g. Ribo 2005) to X-rays (e.g. Corbel et al. 
2002). At higher energies, gamma-ray radiation could also 
be produced in jets (e.g. Bosch- Ramon et al. 2006). A TeV 



Send offprint requests to: Anabella T. Araudo: 
aaraudoQfcaglp. unlp. edu. ar 
* Fellow of CONICET, Argentina 
** Member of CONICET, Argentina 



flare was detected by MAGIC from the HMMQ Cygnus X- 
1 (Albert et al. 2007). Transient gamma- ray events may 
also have been detected from the high-mass X-ray bina- 
ries LS 5039 and LS 1+61 303 by HESS (Aharonian et al. 
2005) and MAGIC (Albert et al. 2006), respectively as sug- 
gested by Paredes (2008). In addition, the HMMQ Cygnus 
X-3 might have been observed flaring in the GeV range 
by AGILE (ATels 1492, 1547 and 1585, see however Atel 
1850). Both this instrument and Fermi also found several 
transient GeV sources in the Galactic plane without known 
counterpart (ATel 1394, Abdo et al. 2009). All this strongly 
variable gamma-ray emission may have a similar origin. For 
instance, the interaction between the jet and the stellar 
wind of the companion star could produce TeV flares (e.g., 
Romero et al. 2003, Romero & Orellana 2005, Romero et 
al. 2007, Albert et al. 2007, Perucho & Bosch-Ramon 2008, 
Owocki et al. 2009). 

In this paper, we propose a model to explain these 
gamma-ray flares, based on the interaction between the 
jets of a HMMQ with wind inhomogeneities. The clumps 
can eventually penetrate the jet, leading to transient non- 
thermal activity that may release a significant fraction of 
the jet kinetic luminosity in the form of synchrotron, inverse 
Compton (IC), and proton-proton (pp) 7r°-decay emission. 

This work is organized as follows. In the next section, 

we describe the main characteristics of the scenario adopted 
and estimate the relevant timescales of the jet-clump inter- 
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Fig. 1. Sketch of a HMMQ with clumpy stellar wind, adapted 
from Romero et al. (2007). 

action; in Sect. 3, we study the acceleration of particles and 
estimate the non-thermal radiative losses; in Sect. 4, we 
describe the calculation of the non-thermal emission and 
present the main results; and finally, in Sect. 5, we draw 
some conclusions. CGS units are used consistently through- 
out the paper. 

2. The physical scenario 

To study the interaction between a clump of the stellar wind 
and a jet in a HMMQ, we adopt a scenario with similar 
characteristics to the binary system Cygnus X-1. We fix 
the separation between the compact object and the massive 
star to be a = 3 X f 0^^ cm (0.2 UA). For the star luminosity 
and temperature, we assume that — 10^^ erg s~^ and 
= 3 X 10'* K, and the stellar mass loss rate is adopted to 
be = 3 X 10~^ Mq yr~^, with a terminal wind velocity 
Uw 2.5 X 10* cm . A sketch of the scenario is presented 
in Fig. 1. 

2.1. Clump model 

The clump is assumed to be spherical and homogeneous. 
Given the uncertainties in the clump parameters, we con- 
sider two values of its size: i?c = 10*" and 10** cm (i.e., 
~ 3 X (10^^ — 10^^) a). For effective jet penetration, a 
large density contrast between the clump and the jet is 
required. Therefore, we assume dense clumps with ric = 
10*^ cm"'', which correspond to a porous clumpy wind with 
a filling factor (clump versus interclump volume ratio) of 
/ = Mi, / Ana'^mpVcnc = 0.005, where nip is the mass of the 
proton* . We assume that the velocity of the clumps equals 
the velocity of the wind, i.e., Vc = 2.5 x 10* cm s~*. The 
temperature of the clumps is taken as Tc = 10* K (Krticka 

* See Owocki & Cohen (2006) for the concept of a porous 
wind. 



Table 1. Adopted parameters in this work. 



Parameter description [unitsj 


Clump 


Jot 


Radius [cm] 


10*^' - 10** 


1.5 X 10'' 


Velocity [cm s~*] 


2.5 X 10* 


10*" 


Density [cm~"^] 


10*2 


4.7 X lO'^ 


Binary system 


System size [cm] 


3 X 10'^ 




Star luminosity [erg s~*] 


1039 




Star temperature [K] 


3 X 10* 




Mass loss rate [M© yr~*] 


3 X 10"'' 




Wind velocity [cm s~*] 


2.5 X 10* 





& Kubat 2001), which is moderately lower than the tem- 
perature at the surface of the massive star. 

2.2. Jet model 

Radio observations demonstrated that the jets of MQs are 
strongly coUimated (Miller- Jones, Fender & Nakar 2006). 
We assume that the jet radius is one tenth of the jet 
height, i.e. Rj{z) — 0.1 z, which corresponds to a semi- 
opening angle oi (f) = 6°. The expansion velocity results 
in Vexp = 0.1 Vj, where Vj is the velocity of the jet. This 
expansion velocity for a free-expanding hidrodynamical jet 
implies a jet base Mach number Vj/v^xp ^ 10, i.e. a 
strongly supersonic outflow. We consider a jet dynamically 
dominated by cold protons with a weakly relativistic bulk 
velocity Vj = 0.3 c, i.e. a Lorentz factor 7j = 1.06. We 
neglect the curvature of the jet produced by its interac- 
tion with the stellar wind. In HMMQ jets of luminosity 
> 10^^ erg s~*, the geometry should not be strongly mod- 
ified (Perucho & Bosch- Ramon 2008), although this effect 
could be important in systems such as Herbig-Haro jets 
interacting with a stellar wind, as studied by Raga et al. 
(2009). 

The kinetic luminosity of the jet is taken to be Lkin ~ 
3 X 10'^^ erg s^*, similar to that of Cygnus X-1 (e.g., Gallo 
et al. 2005, Russell et al. 2007). Using the equation 

-^kin = o"j(7j - l)mp rij Vj, (1) 

where ctj — ttR? is the cross-section of the jet, we can es- 
timate the density nj of the jet material in the laboratory 
reference frame (LRF). At the height of the jet-clump in- 
teraction, which is taken to occur at z = a/2, we obtain 
Tij = 4.7 X 10'' cm"''. Thus, the ratio of the clump to the jet 
densities is x = 2.1 x 10*. This parameter will be very rele- 
vant to the jet-clump interaction estimates. The parameter 
values of the clump and the jet are shown in Table 1. We 
assume an ideal gas equation of state, P = nkT, for both 
the clump and the jet, where n is the gas density. 

2.3. Dynamics of the interaction 

Clumps are formed in the wind acceleration region within 
one stellar radius of the surface of the star (e.g. Puis et al. 
2006) and some of them reach the jet/wind interface. The 
very large inertia of the clumps, linked to the exceptionally 
high density contrast %, allows them to cross the boundary 
of the jet and fully penetrate into it. A large value of x 
is also required to ensure that the clump is not strongly 
affected when penetrating into the jet. 
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To study the physical processes of the interaction, we 
consider the coUision of a single clump with the jet. For 
simplicity, we assume that, on the relevant spatial scales 
of the interaction, the jet is cylindrical. In the context of 
this work, thermal conduction, clump expansion, magnetic 
fields and gravitational forces are not dynamically relevant 
and will be neglected. 

In the LRF, the clump will take a time tc, or jet pene- 
tration time, to fully enter the jet. This provides the first 
relevant timescale 



tr. 



2Rc/vc, 



(2) 



which is tc ~ 80 and ~ 800 s, for = 10^° and 10" cm, 
respectively. In addition, the clump crosses the jet roughly 
at the wind velocity in the jet-crossing time 



2R; 



1.2 X 10^ s. 



(3) 



From the moment when the clump interacts with the jet, 
the ram pressure exerted by the latter produces a shock in 
the clump, which propagates in the direction of the jet mo- 
tion. Assuming that a significant fraction of the momentum 
flux is transferred to the clump, the clump-crossing time, 
which is the characteristic timescale of the jet-clump inter- 
action, can be defined to be 



tr. 



2 Rc 2 i?eVX 



where 



(4) 



(5) 



is the velocity of the shock moving through the clump, 
derived by equating the jet/clump ram pressures. The 
timescale tec is ^ 3 x 10^ and 3 x lO'^ s for i?c — 10"^° 
and 10^^ cm, respectively. 

A shock (the bow shock) is also formed in the jet when 
its material collides with the clump. We assume that the 
bow-shock region width (or clump/bow-shock separation 
distance) is a; ~ 0.2i?c (van Dyke & Gordon 1959), and 
thus the time required to reach the steady state regime is 



i).2Rr 



(6) 



where wj ps ~ Wj/4 is the jet postshock velocity in the LRF. 
Therefore, ths ~ {h/2)tcc/ y/x, i-S- ^bs < ^cc- In Fig. 2, a 
sketch of the situation is shown. 

Once the clump is inside the jet, the latter transfers 
momentum to the clump, accelerating it to the background 
velocity, Vj . The acceleration can be obtained from the force 
exerted by the jet ram pressure and the clump mass (e.g.. 
Fragile et al. 2004): 



9 ■ 



xRc 



(7) 



The jet accelerates the clump to the velocity Vj in a time 



t^ 



(8) 



which is much longer than tec- When one fluid exerts a force 

against another fluid of different density, the hydrodynam- 
ical Rayleigh- Taylor (RT) instability eventually develops. 




Bow-shock 



Jet 



Clump 




Massive 
Star 



Compact Object 
Fig. 2. Sketch of the jet-clump interaction. 



leading to the perturbation and potential disruption of the 
clump. The timescale for this instability is 



i/g, 



(9) 



where I is the instability length in the perturbed region. 

After the bow shock is formed, the jet material sur- 
rounds the clump and we have two fluids that have a large 
relative velocity. This situation leads to Kelvin-Helmholtz 
(KH) instabilities. The timescale for a KH instability to 
develop is 



tKH 



Wx 



(10) 



where Wrd is the relative velocity between the jet material 
and the clump material (we assume Vrei ~ v^). In addition, 
taking into account that we are interested in instability 
lengths similar to the clump size, i.e. I ^ Rc, we find that 
iRT ~ ^KH ~ ice- 

According to the timescalcs estimated in the previous 
paragraphs, the clump can fully enter into the jet if tc ~ 
tcc/5. At this stage, we do not consider the penetration 
of the clump into the jet, and assume that the former is 
completely inside the latter (i.e., the system has cylindrical 
symmetry) . The bow shock is formed in a time much shorter 
than tc and tec {%& ~ tcc/^Jx)- We note that the clump 
might not escape the jet if > tRx/KH- In such a case, the 
clump will be destroyed inside the jet. However, numerical 
simulations show that the instability timescales are longer 
than the clump crossing time by a factor of a few (e.g., 
Klein, McKee & Colella 1994), i.e., tRT/KH > ^cc- 

Regarding the shock properties and given the particu- 
lar characteristics of our scenario, the shock in the clump 
is strong, radiative, and slow, whereas the bow shock is 
strong as well, but adiabatic and fast. For these reasons, 
the shocked and heated material of the clump will radiate 
a non-negligible part of the energy tranferred by the jet. 
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2.3.1. Thermal emission from the clump 

To estimate the density, nps(C), and temperature, Tps{C), 
of the shocked clump at a distance ( from the shock, we as- 
sume pressure equihbrium between the clump and the jet 

shocked plasma and use the equation of energy flux conser- 
vation, taking into account the main channels of thermal 
radiative losses: A(T) 
timescale of thermal cooling equals 



7 X 10-19 erg cm'^ s'K The 



ith = 3 X 10 



2TUO 



1.6 



^ps(C) 



(11) 



The temperature and density in the adiabatic zone of the 
post-shock region are 8.5 x 10^ K and 4 x 10^^ cm"^, re- 
spectively. From these values, the shocked clump material 
cools in a time tth ~ 10 s, which is shorter than tec, which 
means a distance Cth ^ ith^'cs/4 ^ 2 x 10* cm from the 
shock front. Since C,th < 2i?c, the shock in the clump is 
radiative. At distances greater than (th, the density of the 
shocked matter can reach values as high as 10^'' cm~^. 

Although we are mainly interested in the high-energy 
emission produced by the jet-clump interaction, we esti- 
mate for comparison the free-free emission generated by the 
shocked and heated clump. Considering Tpa{(^) and nps(C), 
we estimate the free- free luminosity by integrating the emis- 
sivity (Lang 1999) along the shocked clump, as demon- 
strated in Zhekov & Palla (2007). The corresponding ther- 
mal luminosities are Lth ~ 5 x 10'^° and 5 x 10^^ erg 
for Rc = 10^^ and 10^^ cm, respectively, peaking in the 
soft X-rays. In Fig. 7, the thermal radiation of the clump 
is shown together with the non-thermal emission. 

In contrast to what occurs in the clump, the bow shock 
is adiabatic and fast. For this reason, it is a propitious place 
to accelerate particles to relativistic energies. 



3. Particle acceleration and non-thermal cooling 

3.1. Particle acceleration 

In the presence of a shock and magnetic field, non- 
relativistic diffusive (Fermi I) shock acceleration can occur 
(e.g. Drury 1983). In the linear limit of this theory and 
considering Bohm diffusion, electrons and protons will be 
accelerated to an energy Ef.,p (where e stands for electrons 
and p for protons) in a time 



t. 



E. 



rj 



qBc' 



(12) 



where t] ~ (8/3)(c/ws)^ for perpendicular shocks (Protheroe 
1999). In Eq. (12), Vg is the shock velocity (taken to be equal 
to fj), B is the magnetic field in the acceleration region, c 
is the speed of light, and q is the electron charge. Thus, to 
obtain efficient acceleration of particles, i.e., a short face, a 
high value of B and a strong and fast shock are necessary. 

To study the jet-clump interaction, we consider two val- 
ues of the magnetic field in the jet shocked (bow-shock) re- 
gion, _Bbs- First, we consider a value of the magnetic field 
obtained by assuming that the magnetic energy density is 
10 % of the plasma internal energy density downstream. 
We adopt a sub-equipartition value to ensure that -Bbs is 
dynamically negligible with respect to matter. Fixing 



O.lUr, 



(13) 



where Up is the energy density of the shocked jet gas, given 
by 



(14) 



we obtain Sbs = 150 G. This magnetic field, plus Vs = Vj, 
yield <acc ~ lO~^Ee,p s. 

In addition, we adopt Bbs = 1 G to check the impact 
on our results of a magnetic field significantly weaker than 
in the sub-equipartition case. 

3.2. Non-thermal cooling processes 

Relativistic leptons lose their energy by different non- 
thermal radiative processes, such as synchrotron radiation, 
IC scattering and relativistic bremsstrahlung. On the other 
hand, relativistic protons can also lose energy by means of 
pp interactions. Finally, the shocked plasma can produce 
thermal radiation if the density is high enough, as shown 
in Sect. 2.3.1. 

Relativistic electrons that lose their energy by syn- 
chrotron radiation have the following cooling time (e.g. 
Ginzburg & Syrovatskii, 1964) 



^syn — 



4.1 x 10^ 



s . 



(15) 



Synchrotron radiation is the most efficient radiative process 
in the bow-shock region in the case of _Bbs — 150 G, where 

tsynch ~ 2 X 10~^ E~^ S. 

At the interaction height considered in this work, Zint = 
a/2 = 1.5 X 10^^ cm, the energy density of the photons from 
the star is Wph = 2.4 x 10^ erg cm~'^, the typical photon 
energy being eo ~ 10 eV. For y = eo-Be/(5.1 x 10^ eV)^ > 1, 
i.e. Eg > 2.6 X 10^° eV, the IC interaction occurs in the 
Klein-Nishina (KN) regime. A formula for the IC cooling 
time valid in both a Thompson (Th) and KN regime in a 
photon field with a narrow energy distribution is e.g. Bosch- 
Ramon & Khangulyan (2009) 



he 



6.1 X 10^2 eo (l + 8.3y) (1 + l.Sy^) 



s. (16) 



Mph ln(l + 0.2 y) (1 + O.by + 1.3y^) 

For B ~ 1 G, IC is the dominant loss leptonic channel. 

According to the low particle density of the shocked jet, 
ribs = 4nj = 2 X 10* cm""^, relativistic bremsstrahlung 
losses are negligible in the bow-shock region, the cool- 
ing time for a completely ionized medium being (e.g. 
Blumenthal & Gould, 1970): 



^Brem — 



1.4 X 10 



16 



nZ2 (In 



0.36 



(17) 



where Z is the atomic number, and in our case Z = 1. 

Regarding hadronic emission, 7-rays are produced if rel- 
ativistic protons interact with nuclei through inelastic col- 
lisions. As ill the case of relativistic bremsstrahlung, the 
proton cooling time due to pp interactions depends on the 
density n (e.g. Aharonian & Atoyan, 1996) related 



2 X 10^^ 



(18) 



In the bow-shock region, this process is negligible. 
Otherwise, if relativistic protons accelerated in the bow 

shock penetrate into the clump, pp interactions can become 
an efficient process to generate gamma-rays. 
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10 

log(E, [eV]) 

Fig. 3. Acceleration and radiative loss (synchrotron and IC) 
time for electrons in the bow-shock region. The advection and 
diffusion times are shown for Rc — 10^^ (thin line) and lO'^^ cm 
(thick line). This figure corresponds to the case Bha = 1 G. 



3.3. Maximum energies 

Taking into account energy gains and losses, the maxi- 
mum energy achieved by particles accelerated in the bow 
shock can be easily estimated. Concerning energy losses, 
we consider radiative cooling (described above) and the 
escape of particles. The latter takes into account the ad- 
vection of relativistic particles by downstream bow-shock 
material (tadv ~ Rc/vjps) and the diffusion of particles 
(idiff ^ x'^/2Db, where x = 0.2 i?c and Db = Ee^pc/3qB\^s 
is the diffusion coefficient, which is assumed to be the Bohm 
one). The corresponding timescale is 



l{tdiff , tadv} ■ 



(19) 



Given the small thickness of the bow-shock region, we as- 
sume that particles that escape via diffusion go to the 
clump. 

For electrons, the maximum energy is constrained by the 
escape of particles via diffusion (Bbs = 1 G) and by syn- 
chrotron radiation (Sbs = 150 G), as it is shown in Figs. 
3 and 4. From these figures, we determine that the most 
relevant radiative process in the bow-shock region can be 
IC scattering if _Bbs = 1 G, and synchrotron radiation if 
Bbs = 150 G. On the other hand, the maximum energy for 
protons accelerated in the bow shock is constrained by the 
Hillas criterion (Hillas 1984), i.e. when the proton gyrora- 
dius becomes equal to the size of the acceleration region: 



(20) 



In Table 2, maximum energies for electrons and protons of 
different values of -Bbs and Rc are shown. 

4. Production of gamma-rays and lower energy 
radiation 

We calculated the spectral energy distribution (SED) of 
the emission produced by the most relevant non-thermal 
radiative processes. In the bow-shock region, we considered 
synchrotron and IC radiation, and in the clump, we also 
considered relativistic bremsstrahlung and pp. 
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Fig. 4. Acceleration and radiative loss (synchrotron and IC) 
time for electrons in the bow-shock region. The advection and 
diffusion times are shown for i?c = 10^^° (thin Hne) and 10^^ cm 
(thick line). This figure corresponds to the case Bbs = 150 G. 



Table 2. Maximum energies achieved by acelerated particles in 
the bow shock. The values are in eV units and represent the 
different values of _Rc and B studied in this work. 



-Rc [cm] 


10^" 


10^" 


10" 


10" 


Bbs [G] 


1 


150 


1 


150 


B^"" [eV] 


1.5 X 10'' 


8 X 10'' 


1.5 X 10'^ 


8 X 10" 


B™^" [eV] 


6 X 10" 


9 X 10^^ 


6 X 10^2 


9 X 10^* 



4.1. Distribution of relativistic particles 

We assume an injected population of relativistic particles 
(electrons and protons) in the bow-shock region that follows 
a power-law energy distribution of the form 



Qe,p{Ee,p) — Ke^pE^p exp( — i?e,p/i?™p'*^) 



(21) 



The index F is fixed to 2, typical of linear diffusive shock 
acceleration, and we add an exponential cut-off at high en- 
ergies. The normalization constant Kf, p is determined by 
assuming that 25 % of the jet luminosity incident on the 
bow shock, i.e. (CTc/(Tj)Lkin (where (Tc — T^i?^ ^^'^ clump 
effective cross section ^), is converted into power of the ac- 
celerated particles. We note that the predicted fiuxes scale 
linearly with the adopted non-thermal fraction. 

To estimate the particle energy distribution N^.p, we 
solve the kinetic equation in the one-zone model approxima- 
tion for the bow-shock region (e.g. Ginzburg & Syrovatskii 
1964) 



dN, 



d 



dt 



dE, 



e.p 



(22) 



e,p 



where t is the time and Pe.p = —dEe.p/dt is the energy loss 
rate for particles. 

The relativistic leptons reach the steady state well be- 
fore the shock has crossed the clump. As shown in Figs. 3 



^ We neglect here the region where the bow shock be- 
comes strongly oblicuous, and focus on where this shock is the 
strongest, i.e. right in front of the clump. 
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and 4, the most energetic electrons can diffuse to the clump 
(-Bbs = 1 G) or lose their energy inside the bow-shock region 
by synchrotron and IC radiation {B^s = 150 G). However, 
particles downstream with low energies escape advected in 
the shocked material of the jet before cooling radiatively, 
producing a break in the energy spectrum of particles. By 
equating the advection and synchrotron loss times, we can 
estimate the break energy in the case with i?bs = 150 G, 
which results in = 3 x 10^ and Ei, = 3 x 10^ cV, 
for i?c = 10^" and 10^^ cm, respectively. For the case of 
i?bs = 1 G, the break energy is determined by the advec- 
tion and diffusion times, given i?b = 0.3£^™*''. The electrons 
with energies Eg > E^ can reach the clump and radiate in- 
side it. The energy distribution of these electrons is 



where 



N^,{E,) ~ QUE,) tcooi . 



(23) 



where is the part of Qe that corresponds to 0.3 -E™'^^ < 

E, < E^-^, and ieool = l/(tsyn + ^fc + *Brcm + tion)' The 

ionization cooling term, tion ^ 2 x IQ^^ [Ee/n) s, comes 
from the high densities in the clump. 

As noted in Sect. 3.2, relativistic protons do not suffer 
significant pp losses in the bow-shock region {tdis <C tpp). 
These protons can also reach the clump if they are not 
advected by the shocked material of the jet. By assuming 
that tdiff < iadv, the minimum energy necessary to reach 
the clump is 

= 0.025i;^'*^ , (24) 

and the maximum energies of these protons arc determined 
by the Hillas criterion and shown in Table 2. On the other 
hand, to confine relativistic protons in the clump, the mag- 
netic field must be > 10^ G, much stronger than the ex- 
pected clump magnetic field^, i.e., the protons will cross 
the entire clump in a time ^ Rc/c < tpp. We note that i?c 
should be the shocked clump width if the clump has already 
been affected by the shock. The distribution of relativistic 
protons in the clump is 



Np{Ep) = ^Qp(Ep) 



(25) 



With the steady distributions of relativistic electrons in the 
bow shock, Ne{Ee), and of both electrons and protons in the 

clump, N^{E,,) and Np{Ep), respectively, wc can calculate 
the SEDs of the radiation produced by these particles. 

4.2. Non-thermal emission from the bow shock 

Since the energy density of the synchrotron emission is 
lower than that of the magnetic and stellar fields, syn- 
chrotron self-Compton processes will be neglected in our 
calculations. 

Synchrotron radiation is computed using the standard 
formulae given in Blumcnthal & Gould (1970). Assuming 
that the distribution of relativistic electrons in the bow- 
shock region is isotropic and moves with a non-relativistic 
advection speed, we calculate the synchrotron SED as 



P(£e,e) = 103^Bbs(e/i^c)^'/')exp(-e/Ec)s-^ 



(27) 



is the power function. The magnitude Ec{Ef,) = 5.1 x 
10^^ B E^ erg is the characteristic energy of the produced 

photons. 

For a thermal distribution of target photons, we can 
estimate the IC emission from 

eiic(e) = nph c e / ^ Ar^(£^)^^B£l!£i:^rf^^ , (28) 

where da{xic,eQ,Ee)/dEe is a parametrized differential 
cross section for an isotropic target photon field, valid in 

both Th and KN regimes (Blumcnthal & Gould 1970). This 
magnitude depends on the adimensional parameter 



4en 



1 



(29) 



for m1c^/(AEl) < xic < 1, where eo is the energy of the 
target photons. 

We calculated eLsynch(e) and e_Lic(e) for different values 
of i?bs and i?c) and the results are presented in Figs. 5 
and 6. As seen in these figures, the synchrotron component 
is more luminous than the IC one in the cases of Sbs = 
150 G, reaching bolometric luminosities of Lsynch ~ 10^^ 
and 2x lO'^'"' erg s^^ for i?c = 10^° and 10^^ cm, respectively. 
In contrast, for i?bs = 1 G, the dominant radiative process 
is IC scattering, which leads to bolometric luminosities of 
Lie - 2 X 10^^ erg s'^ and ~ 10^^ erg s'^ for = 10^° 
and 10^^ cm, respectively. The maximum energies achieved 
by photons can be as high as E™i^^ ^ 1 TeV. We note the 
spectral break in the synchrotron and IC emission for the 
case Bbs = 150 G. This feature is produced by the advection 
of particles out from the bow-shock region. In addition, the 
impact of KN losses hardening the electron spectrum can 
also be seen in the spectral shape of the synchrotron and 
IC emission for R^ = lO'^^ cm and Bbs = 1 G. 

We accounted for 7 — 7 absorption in the stellar field 
of the massive star. Since we did not focus on a particular 
binary /observer system geometry, we assumed an isotropic 
target photon field (and hence also neglected angular ef- 
fects in the IC calculations). As seen in Figs. 5 and 6, the 
attenuation in fiux can be of several orders of magnitude at 
energies of hundreds of GeV. In only some cases of specific 
geometries in the 7 — 7 interaction, can the attenuation 
be very much reduced (see e.g. Khangulyan, Aharonian & 
Bosch-Ramon 2008). 

In this work, we have focused on the high-energy emis- 
sion. Radio emission is also produced by electrons that leave 
the bow shock well before they lose most of their energy. 
These particles may radiate in the radio band further down 
the jet and their treatment is beyond the context of this 
study. Because of this, we have not considered the effect of 
synchrotron self-absorption in the radio spectrum, and we 
do not make predictions in this energy range. 



eLsynch(f ) 



N,{E,) P{Ee,e)dEe 



A'" 



Assuming equipartition, magnetic fields 
expected. 



(26) 4.3. Non-thermal emission from the dump 

The most energetic particles accelerated in the bow-shock 
100 G would be region can penetrate the clump. In the case of protons, 
uncooled particles with Ep > 0.025i?™'^'' can diffuse up to 
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emissivity of pions, the specific luminosity of the photons 
produced by 7r°-decay is 



-5-3-1 1 3 5 7 9 11 13 
log(£ [eV]) 

Fig. 5. Computed synchrotron (dotted line) and IC (dashed 
line) SEDs for the emission produced in the bow-shock region 
for the case of Rc — W^'^ cm with Bhs = 1 (thin line) and 150 G 
(thick line) . The production and the 7 — 7 absorbed curves at 
very high energies are shown. 
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Fig. 6. The same as in Fig. 5, but for the case Rc — 10^^ cm. 



the clump and radiate only a part of their energy from 
there. On the other hand, for Bbs = 1 G (< i?cquipartition), 
electrons with Ee > Q.3E^^^ will radiate their energy in 
the clump. 

The relativistic protons that diffuse up to the clump col- 
lide with cold protons there, generating both neutral (tt") 
and charged (tt^) pions that decay to 7-rays and leptons, 
respectively. 

The emissivity of tt*^ is given by 



{Ep)NpiEp)dEp , (30) 



where app{Ep) is the cross section of pp interactions es- 
timated by Kelner et al. (2006). With knowledge of the 



eL, = 2e 



zdEj,- , 



(31) 



- o.i£;^'" 

cay of TT^. 



where E^'" = e + ?r4c''/(4e) and E'^''''' - E'j^''''/6. The 
estimated emission is shown in Fig. 7, along with the con- 
tributions of the other processes. The emission reaches lu- 
minosities of as high as Lpp ~ 10^^ erg {Rc = 10^^ cm; 
Fig. 7). For the adopted values of B^s and E^"""^, the 7 - 7 
absorption reduces the final fluxes and can lead to very 
different spectral shapes. 

Secondary pairs with an effective low energy cut-off at 
will be injected inside the clump by the de- 
These pairs radiate most of their energy inside 
the clump by synchrotron, IC scattering, and relativistic 
brcmsstrahlung. In general, primary electron emission from 
both the bow-shock region and the clump (see the next 
paragraph) is more significant than that of this secondary 
component (see Bosch- Ramon, Aharonian & Paredes 2005, 
and Orellana et al. 2007 for discussion of secondary pairs 
in the context of clouds and jets, respectively). Finally, 
we note that very high-energy neutrinos with luminosities 
~ Lpp would also be generated (e.g., Aharonian et al. 2006; 
Reynoso & Romero 2009). 

To estimate the radiation produced by electrons accel- 
erated in the bow-shock, we assume two values of the mag- 
netic field of the clump: Be = 1 and 100 G, the latter being 
similar to the equipartition clump magnetic field. 

The synchrotron and the IC radiation were calculated 
using the same equations given in Sect. 4.2. Regarding 
the relativistic bremsstrahlung, we used the standard for- 
mulae given in Blumenthal & Gould (1970). As shown in 
Fig. 8, synchrotron radiation is higher than IC in the case 
of Be = 100 G, reaching a luminosity ^ 10'^^ erg 
{Rc = 10" cm) at e ~ 0.1 MeV. For B^ = 1 G, on the 
other hand, the absorbed IC component reaches a similar 
luminosity to that of the synchrotron one, ~ 10'^^ erg s~^. 
Relativistic brcmsstrahlung is negligible in all cases. The 
SEDs for Rc ~ 10^^ cm are similar to those shown in Fig. 
8, but the luminosities are oc CTc and therefore two orders 
of magnitude less. We note that the non-thermal emission 
from the clump is similar to, but of slightly lower intensity 
than the bow shock one, being the spectra harder. We note 
that the radio emission produced by leptons in the clump is 
strongly suppressed by ionization losses and the high low- 
energy cut-off of the injected electrons. For clarity, we have 
not plotted the leptonic emission from the clump together 
with the other components, but they can be compared us- 
ing Figs 7 and 8. 



5. Discussion and summary 

We explored the main physical processes and the nature 
of the radiation produced by the interaction between the 
jet of a HMMQ and a clump of the wind of the compan- 
ion star. The penetration of a clump in the jet produces 
two shocks. One shock quickly reaches a steady state in 
the jet, forming a bow-like shock. The other shock prop- 
agates through the clump driven by pressure equilibrium 
in the jet/clump contact discontinuity. The bow shock is 
adiabatic and fast, and particles can be accelerated to very 
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Fig. 7. SEDs for different values of Bbs and Rc; the curves of both absorbed and unabsorbed (thin hues) IC and pp radiation are 
shown. 



high energies, leptons efRciently emitting synchrotron and 
IC radiation. Otherwise, the shock in the clump is slow 
and radiative. It is a non-efficient accelerator but a ther- 
mal emitter as well as a good target for proton-proton col- 
lisions. If electrons from the bow-shock region enter the 
clump, they could also radiate energy efficiently there via 
synchrotron and IC processes. The SEDs of the mentioned 
radiation components have been computed in the context of 
a microquasar with parameters similar to those of Cygnus 
X-1, showing that in some cases the interaction between 
one clump and the jet can yield significant amounts of ra- 
diation. The entire SEDs for different cases of clump size 
and magnetic field are shown in Figs. 7 and 8. 

As noted above, given the limitations of the one-zone 
approximation, the synchrotron emission has been com- 
puted in the optically thin case, neglecting the impact of 
synchrotron self-absorption. Nevertheless, we note that the 
radio fluxes can be quite high, despite probably being self- 
absorbed right at the bow-shock region. Radio emission 
produced further down the jet by electrons accelerated in 
the bow shock may contribute significantly to the radio 



band, i.e., jet-clump interactions should have associated ra- 
dio flares. 

At X-rays, the emission is produced by synchrotron 
(bow shock and clump) and thermal radiation (clump). In 
the cases with Bbs ~ 1 G, the thermal emission, reaching 
luminosities of Lth ~ 10'^^ erg s~^, is of higher intensity 
than that of synchrotron from the bow-shock region, but 
not from the clump if Rc = 10^^ cm. Synchrotron emission 
from the bow-shock region dominates in the X-rays band 
for i?bs = 150 G and isynch ~ lO'^^ erg s~^. In a source 
such as Cygnus X-1, these levels of X-ray emission would 
be overcome by the accretion disk radiation. Nevertheless, 
in the case of fainter jet sources at X-rays, such as LS 5039 
and LS I -1-61 303 (Bosch-Ramon et al. 2007, Paredes et 
al. 2007), the synchrotron X-rays produced during the jet- 
clump interaction may be detectable, and even the thermal 
component may be detectable on the top of a non-thermal 
component in certain conditions (large clumps of relatively 
low densities). 

Inverse Compton scattering in the bow-shock region 
and the clump produces 7-rays up to very high energies 
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Fig. 8. Non-thermal leptonic emission from the clump, for the 
case _Rc = lO'^^ cm and for the two assumed values of Be, 1 
(bottom) and 100 G (top), is shown. The curves for both ab- 
sorbed (thick lines) and unabsorbed (thin lines) IC radiation 
and relativistic bremsstrahlung are shown. 



ity in HMMQs, which are a new tool for studying the winds 
of massive stars. 

Depending on the wind filling factor (or clump density) 
and the clump size/number -A^c- (Owocki & Cohen 2006), 
the flares produced by interactions of clumps with jets in 
HMMQs could be a sporadic phenomenon, for a low num- 
ber of clumps, or may appear as a modulated steady activ- 
ity, for a high number of clumps (see Owocki et al. 2009) . 
For the parameters adopted here, the former case corre- 
sponds to i?c = 10^^ cm and the latter to i?c = 10^° cm. 
Nevertheless, we note that the jet may be disrupted in those 
cases when too many clumps are simultaneously present 
inside the jet. Assuming that jet disruption takes place for 
CTj < iVc X (jc, for the wind and jet properties adopted in this 
work, the jet could be destroyed if i?c ^ 10"'^° cm. However, 
more detailed calculations of the dynamics of the jet-clump 
interaction are required to clarify this issue. 
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that dominate the radiation output in cases of relatively 
weak magnetic fields (e.g., i?bs = 1 G). In our calculations, 
the highest luminosity achieved is L\q ^ 10^^ erg s~^ for 
i?c = 10"'^"'^ cm, although 7 — 7 absorption can reduce the 
emission substantially above 100 GeV. Proton-proton colli- 
sions in the clump can also produce 7-rays at energies that 
may be as high as ^ lO^'' eV (Sbs = 150 G). The maxi- 
mum luminosity obtained by is nevertheless quite mod- 
est, Lpp ~ 10'^^ erg s~^ for = 10^^ cm, although denser 
and/or bigger clumps, and more powerful jets may yield 
detectable amounts of photons outside the 7 — 7 absorp- 
tion range (0.1-10 TeV). We recall that specific geometries 
of the binary/observer system plus a high-energy emitter 
far from the compact object may render the attenuation of 
7-rays much smaller (e.g., Khangulyan et al. 2008). 

As a consequence of the characteristics of the interac- 
tion, the expected emission is transient. The flare duration 
is related to the permanence of the clump inside the jet 
dependent on RT and KH instabilities, which can destroy 
the clump. Since the clump is accelerated inside the jet, it 
may be disrupted after several dynamical timescales. If the 
clump were not destroyed, it would eventually leave the jet 
as a thin and moderately hot slab of plasma, since the shock 
in the clump is radiative. Given the typical dynamical and 
crossing timescales, the entire event may have a duration 
of between a few minutes and hours. 

The flares of jet-clump interactions would have associ- 
ated lower (synchrotron, thermal emission) and higher en- 
ergy components (IC, pp emission), which should not be 
correlated with the accretion disk activity. The total level of 
emission, the importance of the different components, and 
the duration of the flares, can provide information about 
the jet power and size, clump size and density, and mag- 
netic fields in the interaction regions (Romero et al. 2007). 
Therefore, besides the jet itself, the clump properties can 
be probed by observations at high and very high energies 
(and probably also at radio frequencies) of transient activ- 
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